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Abstract 

Q ' Short-range quark-quark correlations are introduced into the quark-meson cou- 

pling (QMC) model in a simple way. The effect of these correlations on the struc- 



ture of the nucleon in dense nuclear matter is studied. We find that the short-range 



correlations may serve to reduce a serious problem associated with the modified 
quark-meson coupling model (within which the bag constant is allowed to decrease 
with increasing density), namely the tendency for the size of the bound nucleon to 
increase rapidly as the density rises. We also find that, with the addition of cor- 
relations, both QMC and modified QMC are consistent with the phenomenological 
equation of state at high density. 
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About a decade ago, Guichon |IJ proposed a relativistic quark model for nuclear mat- 
ter, where it consists of non-overlapping nucleon bags bound by the self-consistent ex- 
change of scalar (a) and vector (cu) mesons in mean-field approximation (MFA). This 
model has been further developed as the quark-meson coupling (QMC) model, and ap- 
plied to various phenomena in nuclear physics (for recent reviews, see Ref.[[|). Recently, 
Jin and Jennings |§ have proposed an alternative version of QMC (called the modified 
QMC, MQMC), where the bag constant is allowed to decrease as a function of density. 

So far, the use of the QMC model has been limited to the region of small to moderate 
densities, because it has been assumed that the nucleon bags do not overlap. It is therefore 
of great interest to explore ways to extend the model to include short-range quark-quark 
correlations, which may occur when nucleon bags overlap at high density. In this paper we 
will introduce these short-range correlations in a very simple way, and calculate their effect 
on the quark structure of the nucleon in medium. We refer to this model as the quark- 
meson coupling (or modified quark-meson coupling) model with short-range correlations 
(QMCs (or MQMCs)). 

Let us consider uniformly distributed (iso-symmetric) nuclear matter with density 
Pb- At high density the nucleon bags start to overlap with each other, and a quark 
in one nucleon may interact with quarks in other nucleons in the overlapping region. 
Since the interaction between the quarks is short range, it seems reasonable to treat it 
in terms of contact interactions. An additional interaction term of the form, £i nt ~ 
' l Pq('i)F a 4'q(i)4>q(j)F a 'ip q (j), may then be added to the original QMC Lagrangian den- 
sity f2j. Here ip q (i) is a quark field in the i-ih nucleon and T a stands for 1,75,7^,757^ or 
&nv (with or without the isospin and color generators). (For the present we consider only 
u and d quarks.) In MFA most of these terms vanish in a static, spin saturated, uniform 
system because of rotational symmetry, parity etc. We shall retain only the dominant 
MFA contributions, namely the scalar- and vector-type interactions: T a = 1 and 70. 

Next we consider the probability for the nucleon bags to overlap, using a simple geo- 
metrical approach. Let us first consider a collection of rigid balls, with a radius R c . In 
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the close-packed structure of nuclear matter we find that the effective volume per ball 
is given by V c = 4a/2-R^. The corresponding density, p c , is then given by the inverse of 
V c , and hence the radius of the rigid ball is related to the density as R c = l/(4y / 2p c ) 1//3 . 
Returning to our problem, we see that for a given nuclear density p#, if the nucleon 
bag radius, R, which is given by solving the nuclear matter problem self-consistently, is 
larger than R c (= 1/(4a/2pb) 1 ^ 3 ) the nucleons will overlap. If R < R c , there is no overlap. 
Of course, one could build a more sophisticated model, allowing for nucleon motion and 
nucleon-nucleon correlations || , but we believe that the present model is sufficient for an 
initial investigation. 

Now consider two nucleon bags separated by a distance d in nuclear matter. They will 
overlap for d < 2R and the common volume is then given by V ov = Vjv(1— 3y/4+y 3 /16) ||, 
where Vn is the nucleon volume (= 47ri? 3 /3) and y = d/R. It is natural to choose the 
probability of overlap, p, to be proportional to V ov /Vn' 

P(2/)ocl-f + ^. (1) 

Of course, this choice is quite model-dependent. In principle we could use an arbitrary, 
smooth function, which goes to unity at y = and zero beyond y = 2 and which respects 
the three dimensional geometry of this problem. In this exploratory study, we take this 
simple form as an example. 

In mean-field approximation the Dirac equation for a quark field in a nucleon bag is 
given by 

[* 7 • d - (m, " 9> + ~ is/fr + f"M%))l^ q = 0, (2) 

where m q is the bare quark mass, a and uj are the mean-field values of the a and uj 
mesons and g% and are, respectively, the a- and u;-quark coupling constants in the 
usual QMC model 0. The new coupling constants, /S„), have been introduced for the 
scalar (vector)-type short-range correlations, and are given by (see Eq. ([]])) 
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We have also taken y = d/R, with d the average distance between two neighbouring nu- 
cleons at a given nuclear density ps - i.e., as explained above, d = 2R C = 2/(4v / 2ps) 1 ^ 3 - 
Note that since the coupling constants have dimension of (energy) -2 we introduce new, di- 
mensionless coupling constants, /| and f% by dividing by the free nucleon mass (M = 939 
MeV) squared. (If the coupling strength is positive the correlation gives a repulsive force.) 
In Eq.(||), (ipqipq) and (i/jfyq) are, respectively, the average values of the quark scalar den- 
sity and quark density with respect to the nuclear ground state, which are approximately 
given by the values at the center of the nucleon in local density approximation || (we 
will revisit this later). 

Now we can solve the Dirac equation, Eq.fl2|), as in the usual QMC, with the effective 
quark mass 

m\ = m q - g q a a + f q s $ q if> q ), (4) 

instead of the bare quark mass. The Lorentz vector interaction shifts the nucleon energy 
in the medium M: 



e(k) = VM* 2 + P + 3(£> + f*(1>ty q )), (5) 

where M* is the effective nucleon mass, which is given by the usual bag energy 

= ^£ + i Si2 a (6) 
R 3 K ' 

Here B and z are respectively the bag constant and usual parameter which accounts for 
zero-point motion and gluon fluctuations 0. The quark energy, Q (in units of 1/-R), is 



defined by yx 2 + (Rm*) 2 , where x is the lowest eigenvalue of the quark, which is given 
by the usual boundary condition at the bag surface 0. 

The total energy per nucleon at density ps is then expressed as 



E, 



tot 



(27T)3p s 



dhjM" + k 2 + 3(<?> + f q M%)) 



+ ^(m 2 a 2 -m> 2 ), (7) 

where kp is the Fermi momentum, and m a and m w are respectively the a and u meson 
masses. The uo field created by the uniformly distributed nucleons is determined by baryon 



number conservation: u = 3g5P-B/ m S = g^PB/mt (where g u = 3g£), while the a field is 
given by the thermodynamic condition: (dE tot /da) = 0. This gives the self-consistency 
condition (SCC) for the a field §: 

4 fdM*\ fte M* 



{2Tifml \ da J J + p 

where 

/ f>M*\ 

-3g q a S N (v) = -9*C N (a). (9) 



\ da J 

Here g a = 3^5 , at(0) and Cjv(ct) = S'jv(cr) /SV(0), with the quark scalar charge defined by 
Sn{&) — fbagd? 'ifiqipq- We should note that because the scalar-type correlation does not 
directly involve the a field the SCC is not modified by it. However, the correlations do 
affect the a field through the quark wave function. 

In actual calculations, the quark density, {^ q ip q ), in the total energy, Eq.flT]), may be 
replaced by 3p^, and the quark scalar density, contributing to the effective quark mass, 
m*, is approximately given as {^> q ip q ) = {m^/ g a )a because of the SCC (see also Ref. 0). 

Now we present the numerical results. First, we choose m q = 5 MeV and the bag 
radius of the nucleon in free space, R , to be 0.8 fm. We calculate the matter propereties 
using not only QMC but also MQMC. For the latter we take a simple variation of the bag 
constant in the medium to illustrate the role of the short-range correlations: (B / Bq) 1 ^ = 
exp(—g^a/M) with g^ = 2.8 and the bag constant in free space, B . In both models, 
the bag constant (in free space) and z parameter are determined to fit the free nucleon 
mass with R = 0.8 fm. We find B 1 ^ = 170.0 MeV (in QMC, B = B at all densities) 
and z = 3.295. The coupling constants, g a and g^, are determined so as to reproduce the 
binding energy (—15.7 MeV) at the saturation density (p = 0.15 fm -3 ). We find that g 2 a 
= 67.80 and gfr = 66.71 for QMC and g 2 a = 35.69 and g* = 80.68 for MQMC. Note that 
the matter properties at po m both models with short-range correlations are identical to 
those of the original models |7|, |J. This is because, in our simple geometric approach, the 
effect of nucleon overlap starts beyond po (see below). 

In Figs.|l| and we present the total energies per nucleon for QMC and MQMC, 
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Figure 1: Energy per nucleoli for symmetric nuclear matter. The dashed curve is the 
result of the original QMC The solid curve is for the QMC with short-range correlations 
(QMCs). The region enclosed with the dotted curves is the empirical equation of state ||. 

respectively. We determine the coupling constant, so as to reproduce the empirical 
value of the energy around Pb/Po = 2.5 ~ 4 ||. This yields the value = 300 for QMCs 
and f£ = 10 for MQMCs. (Note that in MQMCs the overlap probability is much larger 
than that in QMCs at the same density because the bag radius in MQMC increases very 
rapidly at finite density. This is the reason why the strength of /J for MQMCs is much 
weaker than that in QMCs.) For the strength of the scalar-type correlation, since we have 
no definite guideline, we take the same value in both models: fj = 200 (as an illustration). 
From the figures we can see that the nucleon overlap starts around Pb/po ~ 2.7 or 1.3 
for QMCs and MQMCs, respectively. The empirical energies at high densities (the region 
enclosed with the dotted curves in Figs. 1 and 2 |§) are well reproduced in both models 
(in particular, in MQMCs) if the short-range correlations are considered. 

In FigJ^, we show the change of the nucleon mass in matter. We can see that the effect 
of the short-range correlations on the mass is not strong. Correspondingly, the strength 
of the a field in matter is not much altered by the correlations. 
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Figure 2: Same as in Fig.[l] but for MQMC. The dashed curve is the result of the original 
MQMC. The solid curve is for the MQMC with short-range correlations (MQMCs). (In 
the original MQMC the solution of the quark eigenvalue beyond Pb/po ~ 3 cannot be 
found P|.) 



In FigfJ, we present the variation of the quark eigenvalue. In QMCs, the effect of 
the short-range correlations is weak, while in MQMCs the effect becomes very large as 
the nuclear density increases and the nucleons overlap more and more. Since we chose a 
repulsive scalar-type correlation, the effective quark mass in QMCs or MQMCs becomes 
larger than that in the original model as the density grows. This leads to a larger eigen- 
value in QMCs or MQMCs. As a consequence of this repulsive correlation a solution for 
the quark eigenvalue in MQMCs can be found even beyond Pb/po ~ 3. 

Turning next to the size of the nucleon itself, as measured by the bag radius, we see 
in Fig.[5]that the effect of the short-range correlations can be very significant. While the 
effect is small in QMCs, in MQMCs the bag radius starts to shrink as soon as the nucleons 
begin to overlap. We find a similar effect on the root mean square radius of the quark 
wave function. In the original MQMC it is well known that there is a serious problem 
concerning the bag radius. In particular, it grows rapidly at high density [|3], |j, |10[ because 
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of the decrease of the bag constant. However, as we can see from the figure, the inclusion 
of a repulsive (scalar-type) short-range correlation yields a remarkable improvement for 
the in-medium nucleon size in MQMC. 

In summary, we have studied (in mean-field approximation) the effect of short-range 
quark-quark correlations associated with nucleon overlap. We have found that the em- 
pirical equation of state at high density can be very well reproduced using a repulsive 
vector-type correlation. Furthermore, we have shown that a repulsive scalar-type corre- 
lation can counteract the tendency for the in-medium nucleon size to increase in MQMC. 
This may prove to be a significant improvement because there are fairly strong experimen- 
tal constraints on the possible increase in nucleon size in-medium ||. While our inclusion 
of correlations has been based on quite simple, geometrical considerations, in the future 



we would hope to formulate the problem in a more sophisticated, dynamical way [1 LJ and 
to use it to study the properties of finite nuclei (including hyper nuclei [O ) . 
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Figure 4: Ratio of the (lowest) quark eigenvalue in the matter to that in free space (x 
= 2.052). The dotted curves are the results of the original QMC and MQMC. The solid 
curves are for QMCs and MQMCs. 

for the Promotion of Science. 
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Figure 5: Ratio of the in-medium bag radius of the nucleon to that in free space (i?o = 
0.8 fm). The dotted curves are the results of the original QMC and MQMC. The solid 
curves are for QMCs and MQMCs. 
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